Transport properties of a two Co atoms system: a theoretical approach. 
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A system of two interacting Co atoms has been studied in a recent Scanning Tunneling Microscope experi- 
ment by J. Bork et ai, Nature Physics 7, 901 (2011). This kind of system has been intensively analyzed since 
a quantum critical point was theoretically predicted for an antiferromagnetic interacting dimer. We propose a 
two path structure microscopic model capable of explaining, in all the parameter space, the physics of this ex- 
periment. It consists of two atoms interacting directly and indirectly through a hopping between two electronic 
reservoirs that bypasses the atoms. The many-body system is treated in the finite-U Slave Boson Mean Field 
Approximation. Other results obtained using the Logarithmic Discretization Embedded Cluster Approxima- 
tion are also discussed. The system is characterized as being in the Kondo regime for all the parameter space 
experimentally studied, although the effective exchange interaction between the Co atoms is antiferromagnetic. 



PACS numbers: 73.23.Hk, 72.15.Qm, 73.63.Kv 



Introduction. — The two-impurity Kondo model (TIKM) 
and the closely related two-impurity Anderson model (TIAM) 
are two of the simplest strongly correlated models containing 
the essential physics necessary to study a system that under- 
goes a quantum phase transition (QPT) between a phase char- 
acterized by two independent Kondo-screened magnetic im- 
purities and another where the two impurities are locked into 
a dimer, forming a spin singlet state The parameter that 
drives this phase transition is the ratio I /Tj ( , where / is the 
inter-impurity exchange interaction and is the Kondo tem- 
perature of the individual impurities (assumed to be identical). 
Experiments performed in two-impurity systems have 
been able to carry the system from the Kondo screened phase 
to the antiferromagnetic (AF) regime, but without a precise 
control of the exchange interaction between the two impuri- 
ties that is essential to observe the fluctuations associated to 
the quantum critical point (QCP) corresponding to the QPT 
mentioned above. Furthermore, it is well known that both 
the TIKM HH and the Two Channel Kondo Model (TCKM) 
S are easily driven out of their respective QCPs by inter- 
lead hopping perturbations and for this reason it has been very 
difficult to experimentally observe signals of quantum critical- 
ity. To the best of our knowledge, the only known exception is 
the TCKM studied by Potok et al. H3l- Moreover, recent the- 
oretical works have suggested that these classes of QPT s are 
not observable under current experimental conditions ll I 111 211 . 
Nevertheless, alternative setups have been theoretically pro- 
posed 11311 to try and observe them. 



Recently, a remarkable experiment has been performed 
where a Co atom, positioned at the tip of a Scanning Tunnel- 
ing Microscope (STM), is continuously approached to another 
Co atom adsorbed on an Au(lll) surface 11411 . The position 
of the STM tip was varied with sub-picometer (pm) accuracy 
and hence the ratio I/T 1 ^ could be modified almost continu- 
ously. However, as the tip approaches the surface, the direct 
charge transfer between them increases, thus driving the two- 




FIG. 1: (Color online) Schematics of the model studied in this work. 
a and /? represent the two Co atoms. 



impurity system away from the QCP 115 .El. InRef. H.the 
system stayed in a regime away from the neighborhood of the 
QCP, as a peak in the differential conductance (dC), expected 
to appear at the QPT 171. was never observed. This fact 
was interpreted as resulting from the direct hopping between 
the electrodes. Some features observed in the experiments, 
namely the splitting in dC at small tip/surface separations, 
could be explained by a microscopic model that takes into ac- 
count only the hopping between the electronic reservoirs fl 1411 
and which was solved using the Numerical Renormalization 
Group (NRG) method 11811 . This splitting was interpreted as 
originating from the effective exchange interaction between 
the Co atoms EMI. However, other experimental features, 
such as Fano anti-resonances that first become narrower and 
then transform into a peak, before splitting when z is reduced, 
have not been explained using a microscopic model 12211 . Be- 
sides, it is important to physically characterize the system as 
the distance between impurities is decreased. 

The authors show in this work that the double-Co experi- 
ment described above can be interpreted and understood, for 
all inter-impurity distances studied, by a model that also incor- 
porates, as an essential ingredient, a direct hopping between 



2 



the Co atoms. In particular, the results show to what extent 
the interplay between the direct inter-impurity hopping and 
the indirect one influences the transport properties of the sys- 
tem. 




FIG. 2: (Color online) Differential conductance as a function of 
V/Tk for different values of z (separation between the impurities, 
measured in pm, indicated in the right side). Each curve, calculated 
using SBMFA, has been shifted vertically for clarity. 

Model. — The Hamiltonian is written as a sum of three 
terms, namely 



H — if imp + #hyb + -^leads, 



where 



i—a : /3;a 



u 



(2) 



describes the isolated impurities, ej (where i — a, f3) being 
the energy of each localized impurity, U the on-site Coulomb 
interaction, and a = ± is the spin orientation. The contribu- 
tion 

H hyh = ^ t L c\ lr7 C aa + fflCfi.io.Cffo- + tgpC^Cpa- 
cr 

+ t LR C L la C R s a + H.c. (3) 



describes the hybridization of each impurity with the first site 
of its adjacent metallic lead (modeled as semi-infinite non- 
interacting chains, see below), the hybridization between both 
impurities and the direct tunneling between the left (L) and 
right (R) electron reservoirs (in that order). Finally, 
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describes the L, R leads, represented by two semi-infinite 
chains of non-interacting sites with hopping t between ad- 
jacent sites. The model is depicted in Fig. Q] Most of the 
transport properties were calculated within the finite-U Slave 
Bosons Mean Field Approximation (SBMFA) l23ll . although, 
for the sake of comparison, some results were obtained using 
the Logarithmic Discretization Embedded Cluster Approxi- 
mation (LDEC A) l24ll . Finally, for the sake of simplicity, we 
adopt a symmetric model (i.e., if, =tn = t') 12511 . 

We define as the Kondo temperature for each Co atom 
in the two independent single-impurity Anderson models that 
result from taking tLR = t a p = 0, where is obtained 
through the half-width of the impurity local density of states 
(LDOS) at the Fermi energy. Taking, in units of t, U = 0.8, 
t' = 0.25, and e, = -U/2, we obtain T% = 0.0073. Con- 
sidering t = 1.3 eV for Au, one obtains w 9.5 mcV 
1 111 K, roughly the same Kondo temperature measured in the 
experiments with Co atoms on Au 11411 . will be a refer- 
ence for comparison with experiments. The values of U and 
t' were chosen, within the convergence parameter space of the 
SBMFA method, so that U/T - 14, where T = irt' 2 p(E F ), 
being p{Ep) the leads' DOS at the Fermi energy. This U/T 
value assures that the single impurity system is deep inside the 
Kondo regime at zero temperature. 

Our model includes the hopping t^R between the elec- 
tron reservoirs, providing a channel through which the dots 
couple indirectly and also a weaker direct hopping t a p be- 
tween the Co atoms that, as mentioned above, results to be 
essential to reproduce the experimental results. The assump- 
tion that t a p <C thR is in accordance with the hypothesis 
made in Ref. 11411 . namely, that for a vertical approach be- 
tween the STM tip and surface, the interaction between the 
rf-orbitals of the Co atoms is very weak. Assuming the hop- 



(1) pings to follow tj^ji = Ae lz and t a p = Be Sz , where 



z is 

a parameter representing the inter-impurity distance, then a 
decrease in z results in an exponential increase in t a p and 
tLR (for 7, 6 < 0). After an extensive survey, by vary- 
ing the parameters so as to take into account the constraint 
t a /3 << tLR < t, as discussed above, and the obvious con- 
ditions t a p « and tLR << 1 for the largest z- value used, 
we found that the parameters that best allow the reproduction 
of the experimental results are A = 0.4066, B = 0.0305, 
7 = -0.002534, and S = -0.00973 (compare our Fig. 2 
with Fig. 4(b) in Ref. Jl4ll ). For these parameter values, 
varies from 0.19 to 0.47 and t a/3 from 0.0015 to 0.055, in the 
interval —60 < z < 300 (same as in the experiment, in pm 
units). This parametrization allows us to reproduce qualita- 
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FIG. 3: (Color online) Main panel: Conductance G/Go as a function 
of z. When compared to Fig. 4(a) in Ref. iflHl , we see the same 
overall behavior. Inset: renormalized energy level e a = ep — e as a 
function of gate voltage V g . The plateau at the Fermi energy for both 
z = 200 (continuous (black) line) and z = — 60 [(red) diamonds 
line], indicates that the system stays in the Kondo regime. For larger 
hopping values, t a p = 0.06 and Ilr = 0.13, the plateau starts to 
be suppressed, as shown by the (green) dashed curve for z = —150, 
reflecting that the system enters a crossover regime (see text) 



tively and semiquantitatively the experimental results for dC 
for the whole sequence of z values. 

In the SBMFA, dC is calculated, using the Keldysh 
formalism |f2rjj] . For simplicity, without an applied po- 
tential, we assume electron-hole symmetry. The effect 
of an external bias V is taken into account supposing 
that the system is connected to a left and right reser- 
voirs, which Fermi levels are V/2 and — V/2 respectively. 
In this case, the differential conductance can be written 
as, dC = A7rH'Hm{G L (V/2)}lm{G R (V/2)}\Gl R (V/2)\ 2 , 
where Gl = Gr are the reservoirs' non-interacting Green's 
functions and G\ r (lj) is the many-body propagator from 
left to right reservoirs, under the presence of V . To calcu- 
late this Green function, within the SBMFA formalism, the 
renormalized parameters of the Hamiltonian have to be self- 
consistently calculated incorporating the external bias. How- 
ever, it is known 12711 that for a two-impurity system the 
self-consistent results obtained for the equilibrium situation 
(V = 0) are very similar to the non-equilibrium results (fi- 
nite V), as long as V is smaller than a few times and 
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coming a single peak; (Hi) for higher z values, a Fano anti- 
resonance develops centered at V = 0, with increasing width 
as z increases (see also Fig.|4|. The behavior just described is 
qualitatively and semiquantitatively similar to that observed 
in the experiments by Bork et al. (Fig. 4(b) of Ref. 
However, the range of z values for which our results show a 
single-peak is slightly larger than in the experiments, where 
an anti-resonance line shape persists down to z ps —30, while 
our results show a single-peak feature already at z = 0. 

In Fig. |3] we present SBMFA (dark line) and LDECA (solid 
(red) dots) results for the conductance G/Gq as a function of 
z (where Go = 2e 2 /h is the quantum of conductance), which 
can be compared to Fig. 4(a) in Ref. lITill . Since LDECA 
is exact at the Fermi energy l24ll . the excellent agreement 
between SBMFA and LDECA gives support to the SBMFA 
results shown in this paper. The main difference with the 
experimental results is that our G/Gq values increase more 
smoothly as z decreases. This can be associated to the fact 
that, as expected, the experimental results are very dependent 
upon the distance between the atoms. As stated by Bork et al, 
'mechanical relaxation' effects should be at play when a tran- 
sition from 'tunneling' to 'point contact' occurs, as the tip gets 
closer to the surface. Hence, at this transition the real distance 
z and the associated hopping parameters t a p and t^n are dif- 
ficult to determine. Besides, other matrix elements such as ti, 
or tR, assumed to be constant, may vary at this transition. 

The LDOS of the impurities calculated with SBMFA (not 
shown) shows the widening of the Kondo resonance as z de- 
creases. This is in accordance with previous studies in similar 
systems [f30h . In addition, the splitting in dC is directly corre- 
lated to an equivalent splitting in the impurities' LDOS. Note 
that a split-peak in dC was also obtained by using NRG 1 14], 
with a model where only the hopping between the electronic 
reservoirs was included. However, the overall agreement be- 
tween our theoretical results and the experiments requires the 
inclusion, as a crucial parameter, of the direct inter-impurity 
interaction t a p. In addition, using SBMFA, we have calcu- 
lated the phase difference (not shown) between the two chan- 
nels through which the current can flow, namely the channel 
passing directly through the impurities (via t a p) and the chan- 
nel connecting the electronic reservoirs directly (via t^R). We 
obtained that the phase difference is zero for all values of z, 
indicating that these two channels interfere constructively. 

In the SBMFA, when the system is in the Kondo regime, 
there is a plateau at the Fermi energy (Ep = 0) in the renor- 
malized energy level of each impurity (i a — cp = e), as a 



*<tf < tip/U. Under these conditions, the propa- f(mction of gate vohage y Q As can bg geen ^ ±e imet 



gator G\ r (<jS) is almost independent of the external bias V 
Therefore, we will assume its complete independence from V 
and calculate dC as if the system were in equilibrium [28]. 

Results. — In Fig. [2] we present the SBMFA results for dC 
as a function of V/T^. The dC curves present three funda- 
mental features that should be emphasized: (i) for negative 
values of z, there is a double-peak structure, displaying a split- 
ting that decreases with increasing z; (ii) eventually, still for 
negative values of z, the splitting is totally suppressed, be- 



ef Fig- El for the region of interest, the plateau is perfectly 
defined (compare (black) solid curve for z = 200 with the 
almost identical (red) diamonds curve for z = —60), which 
leads to the conclusion that the system remains in the Kondo 
regime for this range of z. It is important to emphasize that, 
although the system remains in the Kondo regime, the split- 
ting in dC is relatively large for the lowest values of z, as 
shown in Fig.|4j discussed below. For larger values of t a p and 
tLR, the plateau is partially suppressed as shown in the (green) 
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FIG. 4: (Color Online) Splitting as a function of z [(black) solid 
curve]. Inset: Width of the anti-resonance as a function of distance 
in units of T%. The SBMFA results are similar to those obtained in 
Figs. 5(b) and 5(a) in Ref. (14}] [here reproduced schematically as 
solid (red) dots]. 



dashed curve for z = —150. Thus, for large hopping values 
the system enters a crossover region, which is compatible with 
measurements presented in Fig. S.6 of Ref. 11411 . transitioning 
from the Kondo to an inter-atomic singlet state. This indicates 
that the direct and indirect couplings between the impurities 
result in an effective antiferromagnetic spin-spin correlation 
between them, suppressing the Kondo state in each impurity. 

The main panel in Fig. [4] shows the dC-splitting (distance 
between peaks) as a function of z (in the interval —60 < 
z < -40). The (black) solid curve represents the SBMFA 
results and the (red) solid dots are the experimental results 
extracted from Ref. ffl . both in units of their correspond- 
ing single-impurity Kondo temperature Tj ( , which is approx- 
imately 75 K w 6.5 meV in the experiments. The (black) 
solid curve shows that the splitting has an almost linear de- 
pendence with z. In the inset of Fig.@] we show the SBMFA 
[(black) solid curve] and the experimental [(red) solid dots 
curve] results for the width of the dC Fano anti-resonance as a 
function of z, in the interval 300 > z > 0, scaled by Tj ( . We 
fit the SBMFA dC curve with a Fano anti-resonance to extract 
its width. The sharper decrease in the dip width for the exper- 
imental results at z ~ 25 pm can be ascribed to the relaxation 
process that takes place, changing the distance between the 
Co atoms, as mentioned above. We stress the fact that as the 
magnitude of the splitting can change when the experiment is 
repeated, a qualitative and semiquantitative description of the 
experimental results should be satisfactory. 

Single impurity case. — In the case of a single Co impurity, 
the experimental dC shows only one dip, which neither dimin- 
ishes its width nor transforms into a peak as the distance be- 
tween tip and surface is reduced 114, 13111 . In order to study the 
differences between the single- and double-impurity cases, we 
have carried out a z dependence study for the single-impurity 



model. The model is depicted in the inset of Fig. [5] The hop- 
ping between the impurity and one of the reservoirs (tji) and 
the hopping between the two reservoirs (£lh) are varied in the 
same way as in the double-atom case. The SBMFA results 
(main panel of Fig. [5]l show an asymmetric anti-resonance, 
as previously obtained i32ll . We checked that the dip's width 
does not decrease by changing z. Besides, the dip reaches 
its minimum at negative values of bias voltage, as noted in 
Ref. 113211 for an 5 = 1/2 impurity, which is our case. We note 
that Co atoms have 5 = 3/2, but, as mentioned in Ref. 13211 . 
there are only quantitative differences between the 5 = 1/2 
and 5 = 3/2 cases. The value of G/Gq at the Fermi level is 
much smaller than in the two-impurity case. Therefore, these 
results show that the reduction of the anti -resonance width, the 
appearance of a peak and its splitting, as was observed in the 
experiments of Ref. fl4Tl . are a consequence of the presence 
of a second impurity, interacting with the first. 
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FIG. 5: (Color Online) Differential conductance for a single impurity 
as a function of V/T%- The behavior is markedly different from the 
case of two atoms, as in the experiments. 



Conclusions. — Summarizing, the model proposed in this 
work to study the transport properties through two Co atoms 
in series correctly describes the behavior observed in STM 
experiments for all the parameter range. In that respect, the 
inclusion of a direct hopping between the atoms, besides the 
one between the electron reservoirs, proves to be an essen- 
tial ingredient. In addition, we find that the electrons interfere 
constructively along the two possible paths. In our model, 
the direct and indirect couplings between the impurities result 
in an antiferromagnetic spin-spin correlation between them. 
This interaction, although producing a splitting in dC, as in 
the experiments, is not enough to take the system out of the 
Kondo regime for the parameter space studied. 
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